Spin-dependent Rotating Wigner Molecules in Quantum dots 
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The spin-dependent trial wave functions with rotational symmetry are introduced to describe 
rotating Wigner molecular states with spin degree of freedom in four- and five-electron quantum 
dots under magnetic fields. The functions are constructed with unrestricted Hartree-Fock orbits and 
projection technique in long-range interaction limit. They highly overlap with the exact-diagonalized 
ones and give the accurate energies in strong fields. The zero points, i.e. vortices of the functions have 
straightforward relations to the angular momenta of the states. The functions with different total 
spins automatically satisfy the angular momentum transition rules with the increase of magnetic 
fields and explicitly show magnetic couplings and characteristic oscillations with respect to the 
angular momenta. Based on the functions, it is demonstrated that the entanglement entropies of 
electrons depend on the z-component of total spin and rise with the increase of angular momenta. 

PACS numbers: 71.35.Cc, 73.21.-b 



I. INTRODUCTION 

The quantum dot systems, especially the few-electron 
dots with precise control of the particle number, have 
been extensively studied in recent years due to their 
abundant application potential and theoretical meaning. 
Experimentally, these systems are investigated as can- 
didates for future quantum electronics, spintronicsi and 
quantum information dcvices2^i3,. Theoretically, such two- 
dimensional (2D) confined systems with or without exter- 
nal fields have significant differences with the extent ones. 
In strong magnetic fields, the 2D electron gas with short- 
range interaction which exhibits the fractional quantum 
Hall effect has been understood based on the Laughlin 
and composite fermion wave functions^i^ It is also antic- 
ipated that there will be the transition from the quantum 
liquid to the Wigner crystal in the magnetic fields corre- 
sponding to very small filling factorsi^'^'^ In quantum 
dots, it has been demonstrated that similar transition is 
much easier than that in extent systems due to the in- 
trinsic long-range interactionsJ^iii And both the exact- 
diagonalization method (ED) and the analytic theory 
confirm that the crystal states in quantum dots are not 
static but rotating ones with certain angular momenta, 
which are often referred to rotating Wigner molecules 
(RWMs). 

In most investigations on few-electron states in quan- 
tum dots in strong magnetic fields, the spins of elec- 
trons are considered full-polarized and it simplifies the 
analytical theory and computational demands. How- 
ever, in recent experiments, it has been possible to fabri- 
cate and study the quantum dots with negligible Zceman 
splittingFi^iii Then the full considerations of spin degree 
of freedom are needed and pose challenge to theorists. 
Although it has been demonstrated that the spin degree 
of freedom does not change the cystallization process in 
quantum dots, it will still bring new physical phenomena 
to the system. 

In this work, we will introduce the four- and five- 



electron trial wave functions for the spin-dependent 
RWMs. Noticing that the similar four-electron wave 
function has been proposed recently by Jain et ali^ and 
their results are consonant with our four-electron case, 
we will demonstrate that the functions highly agree with 
the ED results and concentrate on the different quantum 
characters between four and five electron quantum dots 
without the Zeeman splitting in strong magnetic fields. 
The remainder of the paper is organized as follows. The 
descriptions of the trial wave functions are introduced in 
Sec. II, the discussion of the accuracy of the functions and 
their quantum characters, such as the vortex numbers, 
magnetic couplings and entanglements between electrons 
are presented in Sec. Ill followed by a summary in Sec. IV. 



II. WAVE FUNCTIONS 

In strong magnetic fields, unrestricted Hartree- 
Fock(UHF) orbits describing localized electrons can be 
approximated by displaced Gaussian functions^ 
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where A = ^/h/m*U , fi = ^/^vf+UJf/4: , Ib = ^jTic/eB. 
z ~ a; — iy is the complex coordinate of an electron. 
Zj ~ Xj + \Yj are the centers of the Gaussians and set 
to the positions of the equilibrium positions of electrons. 
The phase factor ensures the gauge invariance of each lo- 
calized orbit. When wo = or tJc S> ^^o, i-e. A = V^Ib, 
u{z, Zj) can be expanded to the first landau level 



u{z,Zj) = ^Q(Zj)<Pi(z) 



(2) 



where ci{Zj) = (Z* /X)' exp{-\Zj\y2\^)/VV.. 

For a dot with A^i spin-up and N2 spin-down elec- 
trons (A^i + N2 = N), if the electrons form a RWM 
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in a single ring configuration, we can construct 
many-body bases |z| 



7T yl \ 



with z- 

component Sz = (A^i — N2)/2 from the locaHzed UHF 
orbits centered at positions Zi ^ Zn in clockwise, i.e. 
Zj = exp(i27r(l — j)/N). t (J,) represents the spin of the 
electron and each position can be only occupied by one 
particle with a certain spin. 

Taking the four-electron case with 5^ = as an exam- 
ple, the many-particle bases are 
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In the second quantization scheme, the N-electron 
Hamiltonian can be written as 



a - aiak- (4) 



ijkl 



The orbital energies and the direct Coulomb integrals of 
UHF orbits which are same for all bases and can be elim- 
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inated from the diagonal elements of H. If we assume 
that the exchange integrals between the neighbor elec- 
trons and next neighbors are Vijji = vi and VikM = ^^2, 
the Hamiltonian will have a simple form. For four elec- 
trons, it is 
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The eigenstates of {H,S'^,Sz) can be obtained by the 
diagonalization of H. The eigenvalues and eigenstates of 
Eq.® are 
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Using Eq.(l2|), the many-particle bases can be expanded to the first Landau level as 
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These many-particle bases are breaking symmetrical and contain the components of all angular momenta. The 
component of angular momentum L can be obtained by projection operator technique^ 
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For four electrons again, they arc (up to a constant) 

\k)L = 

ll+l2+l3+l4. = \L\ / 4 
0<il<;2,0<i3</4 \j = i 

x(_l)('3+WA|^T^;T^;i^;i) 
where fc = 1, 2, and 
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\3)l = -exp(i6^|L|/4)|l)L 

|4)l = -exp(i2^|L|/4)|l)i 

|5)l = -exp(i2^|L|/4)|2)i 
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due to the 7r/2 rotational symmetry of the positions 
Zi ^ Z4. Then six eigenstates with conserved angu- 
lar momentum can be written as the linear combinations 
of |1)l and \2)l using the relations in Eq. (fTO|l . Since 
the phase differences in Eq. pH)) . for the states with a 
certain total spin, only those with proper L have non- 
vanishing combination coefficients. The allowable L and 
corresponding coefficients for the eigenstates are listed in 
TabU It can be also recognized that and |2)l im- 
ply respectively the ferromagnetic and anti-ferromagnetic 
coupling between electrons. This fact will be impor- 
tant in following discussion of the properties of spin- 
dependent RWMs. 
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For five-electron case with = 0.5, repeat the same 
scheme and again for the two bases \Z\, Z\, Z\, zl^) 
and \Zl, Zl, Zl, Z^, Z^) which are ferromagnetic and fer- 
rimagnetic, the components with angular momentum L 
are respectively 

\zl, zl, zj, Z^, Z^)l= 
o<;i</2<i3:0</4<;5 \i=i / 

X (<5'-^ - <5'0 {5'^ - 5'-) <S3('^+'^) (11) 
\z\, z\, z\, Z^, Z^)l^ 

h+h+h+U+h = \L\ f 5 \ "1/2 

o<;i<i2<i3:0</4<;5 \j=i / 

xS'^'^+'^^llljlllliji) (12) 

where (5 = exp(i|7r). And at this time, there are ten 
eigenstates as the linear combinations of the two bases. 
See Tab im for corresponding coefficients. 

TABLE I: Eigenvalues, angular momenta, total spins, combi- 
nation coefficients and magnetic couplings (MC) of the trial 
wave functions for four-electron spin-dependent RWMs 
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TABLE II: Eigenvalues, angular momenta, total spins, com- 
bination coefficients and magnetic couplings (MC) of the trial 
wave functions for five-electron spin-dependent RWMs 
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Si = Aexp(=Fi|7r), S2 = A exp(=Fi|7r), S3 = A exp(=Fi|7r), 
S4 = Aexp(=Fif7r), A = 0.381966. 



It is worthwhile to point out that the projection lead 
to the fact that not all the eigenstates with restored rota- 
tional symmetry are still orthogonal even if the breaking 
symmetrical ones do. It can be seen in TabU that only 
five of six eigenstates are orthogonal. tJj^ and ipi are the 



same state up to a phase factor. The situation is similar 
in five-electron case as shown in TabHIl 

The scheme for constructing the trial functions with 
5^ 7^ 0(0.5) is same. Because in the following of the 
paper we will mainly discuss the RWMs with 5, = 0(0.5), 
the other wave functions are not shown here and will be 
attached in appendix. 

Having got the eigenstates, we can evaluate their ac- 
tual energies by considering the Hamiltonian of quantum 
dots with parabolic confinement and exact interactions 
as 

The energy expectation value (i/'ji/ [?/')/('/' I V') of the trial 
function with angular momentum L are 

EdotiB) = ^huj, + {L + N)h^ + ^{V) (14) 
2 ujc Lb 

where (V) is the expectation value of interaction energy 
between electrons evaluated with Ib equals to a certain 
length ^0- 

Noticing that the overlap of a single-particle orbit ipi (z) 
in a quantum dot with the first Landau level ipi{z)^g=o is 
[2V2AZb/(A2-|-2/^)]I'I+-'^, we can also evaluate the overlap 
between the trial function with angular momentum L and 
the exact-diagonalized one as 

where (Ct, Ce) represents the inner product between the 
determinant coefficient vectors of trial function and that 
of exact-diagonalized one. 

In order to study the quantum correlations in the trial 
functions, we can calculate the von Neumann entropiesi^ 
between an electron and the other part of the system as 

5 = -tr(p/log2pO (16) 

where = {tp\a'^a^\tp) is the single-particle reduced 
density matrix. 

III. DISSCUSSION 

A. Energies, overlaps and vortices of the functions 

The energy expectation values of the trial functions 
obtained from Eq. p^ as function of magnetic field are 
plotted in Fig[TJ In the same plot, we also present the 
energies of lowest states with different total spins calcu- 
lated by ED. The ED calculations are executed beyond 
the lowest Landau level approximation and the confine- 
ment strength is set to 2meV. Then we employ the results 
of ED as a criterion to examine the accuracy of the trial 



4 



> 

CD 



LU 
LU 




B(T) 

FIG. 1: (Color online) Energy expectation values of the trial 
wave functions (solid lines) and the energy of the lowest states 
calculated by ED (dashed lines) for four (a) and five electrons 
(b) with total spins from the minimum to maximum values 
indicated by black, red and green lines, respectively. N times 
of the energy of the first Landau level have been subtracted 
from the total energy. 



functions. It can be seen in the plot that the results 
for four electrons obtained from the trial wave functions 
agree well with the ED ones in the magnetic field stronger 
than 5.7T with only an increase of the total energies lower 
than 0.25meV (0.6% of the total energy at B = 5.7T). 
For five electrons, the difference is lower than 0.398meV 
(still 0.6% of the total energy) when the field is stronger 
than 6T. 

With the change of the field, both trial functions and 
ED exhibit accordant angular momentum transitions for 
the states with different total spins. In fact the allow- 
able angular momenta of the trial functions listed in 
Tab|T] and |TT] are just the angular momenta existing in 
the transitions identified by ED. When the magnetic field 
is strong enough, the trial functions can give correct en- 
ergy sequence of different spin states and the values of the 
field where the transitions take place. In smaller fields, 
the full-polarized rotating Wigner molecules always have 
lower energies than the ones with other total spins if the 
states in quantum dots are described by the trial func- 
tions. However, the ED demonstrates that the states in 
small fields are actually liquidlikc and the states which 
are not full-polarized can have lower energies. Then the 
trial functions cannot describe them correctly. It should 
be noticed that because the field ranges for liquid-to- 
crystal transition depend on the confinements of quan- 
tum dots, the trial functions can be expected to give 
better agreement with ED ones in smaller fields for the 
dot with weaker confinement. 
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FIG. 2: (Color online) Overlaps between the trial wave func- 
tions and ED ones for four (a) and five electrons (b) as func- 
tion of angular momentum, x, • and ■ correspond to the 
states with non-, partial- to full-polarized spins, respectively. 



Besides the energy expectation values, we also eval- 
uate the overlaps between the trial functions and the 
ED ones. As shown in FiglH the higher overlaps con- 
vinced us that the trial functions can describe the states 
with larger angular momenta in quantum dots, i.e. the 
lowest states in strong field. The zero points, i.e. vor- 
tices of the many-body wave functions can reflect the 
charge distributions^ so we also inspect the vortex 
structures of the trial functions. The vortex number 
of the function just equals to the highest angular mo- 
mentum occupied by the single particle. For four elec- 
tron non-full-polarized states with Sz — 0, the vortex 
number generally equals to |L| — 1 as the function con- 
tains the component corresponds to the angular mo- 
mentum distribution {li^hjhjh} = {0,1,0, |X| — 1} or 
{li,l2,h, h} = {0, \L\ — 1, 0, 1}. An exception is the case 
when |L| — 1 = 4n, where n is an arbitrary integer. From 
Eq.([n]) it can be found that the coefficient of such com- 
ponent is zero, so the vortex numbers of these states are 
— 2. For the full-polarized states, the vortex numbers 
are \L\ — 3 due to the Pauli exclusion principle. We also 
compare the vortex structures of the trial functions with 
that of ED ones. Although the angular momentum com- 
ponents of the trial functions are restricted in the first 
Landau level, the vortex number in the scope which con- 
tains the majority of the electron charge can agree well 
with the ED ones. Of course, due to the fact that our trial 
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FIG. 3: (Color online) Spin-dependent conditional probability densities of four-electron states with (L, 5") — (—26, 0) (—27, 1) 
and (-28,1) and five-electron states (—46,0.5) and (—46, 1.5) in QDs.The panels in the first(second) row show the distributions 
of the spin-down(up) electron. The red dots indicate the positions of the fixed spin-up electron. 



functions mainly reflect the long-range limit characters of 
the electrons in quantum dots, the vortex distributions 
are more dispersed and different from the ED ones. 

For the energies, overlaps and vortices, the magnetic 
field or angular momentum beyond which the trial func- 
tions can get agreement with the ED ones for five elec- 
trons is larger than that for four-electron case. This is 
because that the transition from liquid to crystal states 
needs larger fields for the dots with larger particle num- 
bers. 



B. Spin correlations of RWMs 

It has been pointed out in Sec. II that the many-body 
bases and \2)l for four-electron imply two species 
of magnetic coupling existing among electrons, respec- 
tively. Then the spin-dependent RWMs may also have 
corresponding properties since they are linear combina- 
tions of the two bases. 

As listed in TabUl we have indicated the ferromag- 
netic, anti-ferromagnetic couplings and no specific cou- 
pling with number -1, 1 and respectively. There are 
four-electron states like "02^5 which only contain the com- 
ponent or \2)l and of course they should have the 
corresponding type of magnetic coupling. There are also 
states like ipi and tp^ containing both of the components. 
The full-polarized states tpe should not have any specific 
magnetic coupling because the spin parts of the wave 
functions can be separated from the spatial parts. Then 
the states tpi should have residual anti-ferromagnetic 
coupling because the ratio of the two components in them 
is different from that in -06 • 

We also calculate the conditional probability densities 
(CPDs) of finding electrons with a certain species of spin 
to explicitly identify the two kind of magnetic couplings. 
As shown in FigEl if a spin-up electron is fixed, the 
remainder spin-up and -down electrons of four-electron 
state with L —27, 5 = 1 will be clearly at the neigh- 
bor and opposite positions, respectively. There are in- 
deed the probability to find spin-down electrons at the 
neighbor positions. This is because the spin-down elec- 



trons have equal probability to stay at two remainder po- 
sition if the spin-up electron is at left or right neighbor of 
the fixed ones. Such CPDs just reflect the ferromagnetic 
coupling between electrons. Then the other CPDs as the 
case of the state with L = — 28, 5 = 1 are the reflection 
of anti-ferromagnetic coupling. In these two examples, 
the states only contain the single component or \2)l 
respectively. The case of the state with i = — 26, 5* = 
is also anti-ferromagnetic as the states with L = —28, 
5 = 1. However their CPDs have a slight difference since 
the state with L = —26, 5 = is just a sample of ipi 
which contain both components |1)l and \2)l. 

Similar to the four-electron case, the five-electron bases 
Eg. pT]) and Eq. p^ are ferromagnetic and ferrimag- 
netic respectively. Then the five-electron spin-dependent 
RWMs also have corresponding magnetic coupling as 
listed in TabHIl where —1,1 and respectively repre- 
sent ferromagnetic, ferrimagnetic and no specific cou- 
plings. Different from the four-electron case, all of the 
five-electron RWMs contain both components of the two 
bases, they can only have the residual magnetic coupling 
if the ratios of the two components are different from that 
in full-polarized states ^lo- The states like ■01 which have 
same ratio as ipiQ also have no specific magnetic coupling. 
In FigOwe also show the CPDs of two five-electron states 
with L = —46,5 = 0.5 and 5 = 1.5 which respectively 
reflect ferrimagnetic and ferromagnetic couplings. 

Before the finish of the subsection, it is worthwhile to 
point out that there are regular magnetic coupling oscil- 
lations for spin-dependent RWMs with respect to the an- 
gular momentum which are in accordance with the elec- 
tron molecules theoryii^ It can be identified in Tab U and 
|TT]that with the increase of the angular momentum, the 
magnetic couplings of different spin states have periodic 
variations. We argue that these L-dependent spin cor- 
relations originate from that the configuration of local- 
ized UHF centers i.e. the classical position of electrons 
satisfies 2tt/N rotational symmetry. The gauge invari- 
ance of the UHF orbits ensures the phase differences in 
Eg. pT]]) and results in the magnetic coupling oscillations. 
Therefore, we can conclude that these i-dependent spin 
correlations reflect the nature of spin-dependent RWMs. 
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C. Sz- and L-dependent entanglements 
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FIG. 4; (Color online) Entanglement entropies of the trial 
wave functions of four (a) and five (b) electrons as function 
of angular momentum, x , • and ■ correspond to the states 
with non-, partial to full-polarized spins, respectively. 

The rotating Wigner molecules are the states driven by 
strong particle correlation. When spin degree of freedom 
is considered, the correlation is naturally spin-dependent. 
The magnetic couplings discussed in the previous subsec- 
tion arc one kind of the reflection of such spin-dependent 
correlations. 

Now we explore the entanglements of the trial func- 
tions which also reflect the quantum correlations exist- 
ing in RWMs. The von Neumann entropy is used for 
evaluating the entanglement between one electron and 
the other part of the system. In the case of identical 
particles, the exchange symmetry of wave functions con- 
tributes a constant logj N to the total entropies. Such 
parts of the entropies do not reflect the entanglements 
we concern, so wc extract them from the total entropies 
and the corresponding results are shown in Fig[4l which 
not only contain the states with Sz = 0(0.5) but also the 
ones with Sz = 1(1-5) and 2(2.5). The entropies given by 
the trial functions show that the entanglements of RWMs 
increase monotonously with the increase of the angular 
momentum. For a given L, the states with different S 
but same Sz have almost equal entanglement, especially 
when \L\ is large enough. For the states with same S 



but different Sz, the entropies will have a constant dif- 
ference. For the states with largest Sz , since all electrons 
have same spin, there are no entanglements originating 
from spin components and the entropies will differ from 
that of the states with smallest Sz by one. 



IV. SUMMARY 

In summary, a set of trial wave functions for describ- 
ing the spin-dependent rotating Wigner molecular states 
is constructed from the localized unrestricted Hartree- 
Fock orbits with the assumed long-range interaction. The 
projection operator technique is employed to restore the 
rotational symmetry. By examining the energies of the 
functions and their overlaps with the exact-diagonalized 
ones with the full consideration of the interactions, it 
is demonstrated that the trial functions arc suitable for 
the studies of the few-electron states of quantum dots 
in strong magnetic fields. The vortex numbers of the 
functions are generally jLj — 1 for the non- full-polarized 
states with Sz = 0, where \L\ is the angular momen- 
tum of the state, except the states with |L| — 1 = An 
whose vortex numbers are \L\ — 2. They automatically 
satisfy correct angular momentum transition rules of dif- 
ferent spin states with the increase of the magnetic fields. 
These functions explicitly show that there are ferromag- 
netic and anti-ferromagnetic couplings between electrons 
in four-electron rotating Wigner moleculars, and fer- 
romagnetic and ferrimagnetic couplings in five-electron 
case. The different spin states have specific magnetic 
coupling oscillation patterns with the change of angular 
momentum. The allowable angular momenta and the L- 
depcndent magnetic coupling of different spin states are 
caused by the 2tt/N rotational symmetry of the classi- 
cal positions of electrons and reflect the nature of spin- 
dependent rotating Wigner moleculars. It is also shown 
that the entanglements i.e. the quantum correlations of 
the states increase monotonously with the increase of an- 
gular momentum. The states with different S but same 
Sz have similar entanglement entropies. The entropy dif- 
ferences of the states with different Sz are approximately 
constant in strong magnetic flelds. These trial functions 
will be useful in future studies of the few-electron states 
and their quantum behaviors in quantum dots in strong 
magnetic fields. 
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APPENDIX: SITUATIONS FOR / 0(0.5) 

For four-electron RWMs with Sz = I, the many-body 
bases are 




f 1^1 ) ^h^i) \ 

) -^3' -^4' ^2/ 
V 1^2 1 ^3 ' ^4' ^l) / 



(17) 



Examining the CPDs or noticing that there is only one 
spin-down electron, it can be found that the electrons in 
these bases have no specific magnetic coupling. 
The corresponding Hamiltonian is 



H 



-2vi — V2 Vi —V2 Vi 

Vi —2vi — V2 Vi —V2 
— V2 Vi ^2vi — V2 Vi 

Vi -V2 Vi ~2vi ~ V2 , 



(18) 



The eigenvalues and eigenstates of (H, S"^, Sz) are 



5*= 1 



5* = 2 





= -2vi 




■■ (-1,0,1,0) 


E2' 


= -2vi 


V'2' 


: (0,-1,0,1) 


^3' 


= -2V2 


V'3' 


: (1,1,1,1) 


E4' 


= -Avi - 


- 2V2 1p4' 


: (1,-1,1,-1 



(19) 



It can be found that the states with Sz = I have same 
energies as the ones with Sz = and same S. 

Again, we expand |1)' to the first Landau level and 
then the component with angular momentum L is 



|1)'l = 



ll+l2+l3 + l4 = \L\ 



-1/2 
■I 



0<h<l2<hfi<h \ j=l 



x(i'3_i'^)(-i/^|zT,4,4,;i), 

and there are relations as follows 

|2)'^ = exp(i67r|L|/4)|l)i^ 
|3)'^ ^ exp(i47r|L|/4)|l)i^ 
|4)'^ = exp(i27r|L|/4)|l)i^ 



(20) 



(21) 



Using Ea. ((2T|) . four eigenstates in Ea. ((T9|) with re- 
stored rotational symmetry can be obtained. For the 
states with S = 1, only those with \L\ = 4fc±l or |L| = 4fc 
are allowable. For the states with 5* = 2, only those with 
\L\ = Ak + 2 are nonvanishing. These angular momen- 
tum rules are same as the ones for the states with Sz = 0. 
Of course these states do not have any specific magnetic 
coupling. 

Similarly, it can be obtained that the five-electron 
RWMs with Sz = 1.5 have same energies and angular 
momentum rules as the ones with Sz = 0.5. They can 
be represented by the many-body basis (up to a phase 
factor) 



I ^1 7 ^2 ' ^3 7 ^4 ' ^5 ) i — 



-1/2 



li+l2+h+l4+h = \L\ / 5 
0<il<;2</3<i4,0<i5 \i = i / i<j 

XS^''^ l\, l\, Zg). 



(22) 



And they also have no specific magnetic coupling. 

Finally, the wave functions for full-polarized RWMs 
with largest Sz can be found in other references^i and 
are no longer discussed here. 
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